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I Abstract. We report optical, near-infrared and X-ray observations of the afterglow of GRB 0508 14, 

which was seen to exhibit very red optical colours. By modelling its spectral energy distribution we 
find that z = 5.3 ± 0.3. We next present a carefully selected sample of 19 Swift GRBs, intended to 
, estimate in an unbiased way the GRB redshift distribution, including the mean redshift (zmean) as 

■ well as constraints on the fraction of high-redshift bursts. We find that Zmean = 2.7 and that at least 

\^ I 5% of the GRBs originate at z > 5. The redshift distribution of the sample is qualitatively consistent 

. with models where the GRB rate is proportional to the star formation rate in the Universe. The high 

mean redshift of this GRB sample and the wide redshift range clearly demonstrates the suitability 
of GRBs as efficient probes of galaxies and the intergalactic medium over a significant fraction of 
the history of the Universe. 
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. INTRODUCTION 

X 

o3 ] The immense luminosities of the gamma-ray bursts (GRBs), coupled with their origin 

in the core collapse of massive stars IH 0] and their y-ray penetration through dust, 
open up a variety of intriguing cosmological applications. Much effort has been directed 
into the use of GRBs for studying star formation (e.g. [SJ), as backlights for exploring 
high-redshift galaxies and the intergalactic medium (e.g. 0, S), and even as probes 
of cosmological parameters (e.g. ||a,0])- Although the GRB population observed until 
the end of 2004 had enabled much progress in the field, it was widely expected that 
the launch of Swift, and the subsequent order of magnitude increase in the number of 
GRBs open to detailed study, would allow further insight into the high-redshift Universe 
lH]. Indeed, the ability of Swift to locate and follow-up a fainter burst population than 
was previously possible has allowed the study of more distant bursts. The mean 
redshift of pre-Swift bursts was Zmean = 1-4, while we show here that bursts discovered 
by Swift now have Zmean = 2.7, inclu ding the first burst to have been discovered with 
z>6, GRB 050904 at z = 6.295 (e.g. JI^). 
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FIGURE 1. The spectral energy distribution of the GRB 050814 afterglow at Af = 14 hr The strong 
break blueward of the /-band is too strong to be readily explained by reddening alone and is best fit by 
the presence of the Lya break at z = 5.3. The solid curve is a fit to the data at that redshift. The dashed 
line shows the spectral slope expected from a synchrotron emission in the fireball model with j3 = 1 . The 
filter response functions are also shown. The horizontal error bars represent the FWHM of each filter The 
V-band upper limit is 2c7. The inset shows the VRIJK observations (filled squares) along with the X-ray 
spectrum (filled triangles) at At = 14 hr The dashed line is the same j3 = 1 slope as in the main panel. 



SPECTRAL ENERGY DISTRIBUTION OF THE AFTERGLOW 

Our multiband observations of the GRB 050814 afterglow are presented in lITlll . 
They allowed the construction of its spectral energy distribution (SED), displayed in 
Fig. [U where we have corrected the observed data points for foreground (Galactic) 
extinction. The SED has a strong break blueward of the /-band, exhibiting colours of 
I — K = 3.44 ±0.29 mag and R — I = 2.87 ±0.10 mag, corresponding to spectral slopes 
of = 1.78 ±0.12 and pR, = 11. 70 ±0.04, respectively (F oc y-^). The latter value 
is unreasonable for GRB afterglows, implying an electron energy power-law index more 
than ten times higher than normally observed. Even in the case of high local extinction 
(Ay), such steep slopes cannot be obtained (see also [12]). 

The most likely explanation for the steep break observed is due to the presence of 
the Lya break at a redshift of 5 < z < 6. To provide a more robust estimate of the 
GRB 050814 redshift we fit the available photometry at different redshifts, allowing for 
a range in /3 and Ay modelled using the parametrization of [13]. The models of [ 14ll 
provide the average hydrogen opacity as a function of redshift. The minimum is 
obtained for z = 5.3 ±0.3. However, we are only able to obtain weak constraints on j8 
and Ay. Fixing /3 = 1.0, a typical value for GRB afterglows, results in a best fit of a 
restframe Ay = 0.9 mag and an unchanged redshift. This Ay is marginally higher than 




TABLE 1. An update to the list of 28 long-duration GRBs 
from table 2 in lllll ; these additional nine bursts were detected 
after 30 September 2005. Here 0sun is the Sun-to-field distance, 
©Moon the Moon-to-field distance and /Moon the Moon illumina- 
tion at the time the burst occurred. For a burst detected in the 
optical but without a reported redshift, an upper redshift limit is 
estimated based on the filter it is detected in. 
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has been inferred from the SEDs of pre-5w/// bursts with bright optical afterglows 
(OAs), but is a necessary consequence of the red I — K colour. The extrapolated /3 = 1 .0 
line, normalized for Ay = 0.9 mag, slightly overestimates the predicted X-ray flux (inset 
of Fig. [T]), indicating that /3 is a bit steeper; /3 = 1 . 1 would make the X-ray data fall on the 
line. Since the best fit X-ray spectral index is consistent with the assumed optical/NIR 
one, a cooling break between the optical and X-rays can be ruled out. 



THE REDSHIFT DISTRIBUTION OF SWIFT BURSTS 

In order to study the redshift distribution of GRBs, it is important to carefully select 
a sample containing bursts which have "observing conditions" favorable for redshift 
determination. In [11] we introduced four criteria with the aim of constructing such a 
sample. Here we recap those criteria and add a fifth one: (1) Small error circles, hence 
the bursts have to be localised with the XRT. (2) The Galactic extinction in the direction 
to the burst has to be sufficiently small or Ay^^ < 0.5 mag. (3) The XRT error circle 
should be distributed quickly (within 12 hours) for a relatively rapid follow-up. Although 
the automatic slewing of Swift was enabled in the middle of January 2005, part of the 
following month was dedicated to calibration which could not be interrupted. Therefore, 
we have only included bursts occurring after 1 IVIarch 2005. (4) Rejection of bursts with 
a declination unsuitable (above +70° or below —70°) for follow-up observations. (5) 
The Sun-to-field distance has to be large enough, with 0sun 55°. 



The first 28 bursts in the sample were listed in table 2 in lllin . Nine additional bursts 
are presented in Table [B For each burst we have also listed the ]VIoon-to-field distance 
(0Moon) and the IVIoon illumination at the time of the burst. This is done to examine if 
these parameters affect the redshift determination significantly, e.g. a full JVIoon close to 
a burst location. This is of course difficult to quantify as the OA brightness also plays 
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FIGURE 2. The cumulative fraction of GRBs as a function of redshift for 43 pre-Swift bursts (upper 
stepwise curve) and 19 Swift bursts (lower stepwise curve). Oyerplotted are three simple models for the 
expectation of the redshift distribution of GRBs: model II from |21] in which the GRB rate is proportional 
to the star formation rate (solid curve), model IV from [211 in which the GRB rate increases with 
decreasing metallicity (dashed curve) and a model from InM in which the GRB rate is proportional to 
the star formation rate (dash-dotted curve). All three models fold in the Swift/BAJi flux sensitivity. 

a role. For example, GRB 050820A (table 2 in flT']) has a measured redshift even if it 
occurred during a full Moon and 0Moon = 34°. Therefore, we decided not to limit the 
sample further. This relatively "clean" sample of 37 bursts has a redshift recovery rate 
of roughly 50% (19/37). 

Figure|2lshows the redshift distribution of the 19 bursts with a reported redshift in our 
Swift sample. Both the mean and the median is 2.65, more than twice as large as the 
corresponding numbers for pre-Swift bursts. A natural explanation for this increase is 
the lower trigger threshold of Swift compared to previous missions, giving rise to fainter 
(Swift events are on average 1.7 mag fainter in the 7? -band at a similar epoch: 0]) and 
higher redshift bursts. This is complemented by the accurate positions provided by Swift 
and the rapid response of a variety of telescopes aimed at redshift determinations. 

This Swift sample is the most uniform to date and it is of interest to compare its 
redshift distribution to models predicting the fraction of GRBs expected to occur at a 
given redshift. [21] have modelled the expected redshift distribution for GRBs, utilising 
several models including those which follow the globally averaged star formation rate 
(model II), and those which scale according to the average metallicity of the Universe 
at a given redshift (model IV, see e.g. [22, |22|]). ll24l] have also carried out a similar 
exercise, where the GRB rate is assumed to be proportional to the star formation rate. 
These models are plotted in Fig. El 

It is remarkable how similar the observed Swift redshift distribution is to the model 



predictions; we can now reason that GRBs indeed trace star formation (see also ['25', "2^). 
However, with the available sample and the limited flux sensitivity of the 5w«///BAT for 
z > 5 bursts, it is currently not possible to determine if GRBs are unbiased tracers of 
star formation. For example, models II and IV from lETi] are nearly indistinguishable 
when comparing to the relatively small sample of 19 bursts. Note that although model 
II from [21] and the [24] model both presuppose that the ORB rate is proportional to 
the star formation rate, they use different assumptions regarding the poorly determined 
GRB luminosity function and the intrinsic spectral shape, explaining their difference in 
Fig. El 

By including all the bursts in Table [Hand table 2 in [lV\, we are able to constrain 
the number of bursts above a specific redshift. For example, 5%-40% of the bursts are 
located at z > 5. The lEHl and lE^l predictions are 10% and 2%, respectively, suggesting 
that the GRB luminosity function parameters and/or the GRB spectral index assumed 
in 112 ih might be more appropriate. [27.1 also predict that 10% of the Swift GRBs should 
originate at z > 5. 



DISCUSSION & CONCLUSIONS 

The mean redshift of our relatively unbiased Swift sample (Zmean = 2.7) is larger than the 
median redshift of sub-mm galaxies (Zmedian = 2.2: [28]) and is similar to that of Type 
2 AGNs (zmean ~ 3: [29]). With two z > 5 GRBs discovered within a space of a month, 
and primarily due to the spectroscopic redshift of z = 6.295 for GRB 050904 ifioll . we 
are finally accessing the GRB high-redshift regime. Are we starting to probe the era of 
Pop III stars? If the transition between the dark ages and the era of reionization occurred 
around z ~ 6-7 (see e.g. [30] for a review), the answer might be positive. However, [i3lll 
have calculated that at most one massive metal-free star forms per pre-galactic halo, and 
since the GRB progenitor may need to be a member of a close binary system in the 
coUapsar scenario (e.g. ll32Ll33ll34ll ). it seems unlikely that Pop III stars could end their 
lives as GRBs. [35] have proposed a non-binary possibility in the coUapsar scenario, 
introducing unusually rapidly rotating massive stars. It is therefore possible that Pop III 
stars are GRB progenitors, although the number of unknowns is currently too large to 
arrive at a concrete conclusion. 

The sample contains GRB 050814, whose OA was particularly faint in the i?-band; 
the observed optical to X-ray spectral slope is flatter (/3ox = 0.36) than expected for 
the fireball model. Hence, GRB 050814 is classified as a dark burst as defined by ll36ll . 
We have argued that this is most likely due to the high-redshift nature (z = 5.3) of 
this burst; the R — I colour is extremely red which is impossible to explain by strong 
extinction given the observed I — K colour. Indeed, a similar conclusion was proposed 



for GRB 980329 MM 



It is clear that GRBs have now opened up a window to the very high-redshift Universe. 
The emerging GRB redshift histogram (Fig.|2l) strongly indicates that GRBs can be used 
to trace the star formation in the Universe over a wide redshift range (0 < z ^ 7). Future 
instrumentation, such as the X-shooter [38], will hopefully shed light on the end of the 
dark ages and the possible GRB/Pop III connection. 
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